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I. INTRODUCTION 



Magnetic monopoles have raised a large interest since their first construction by Dirac [|TJ . 
While the Dirac-monopole has a singularity, the so-called Dirac string, 't Hooft and Polyakov 
|| came up with the construction of a particle-like magnetic monopole in SU(2) Yang-Mills- 
Higgs (YMH) theory with triplet Higgs scalar. The magnetic charge of this object results 
from the topological properties of the solution and is directly proportional to the degree of 
the map from space-time infinite to the vacuum manifold of the theory. Minimal coupling 
of the SU(2) YMH model to gravity leads (for suitable choice of the boundary conditions) 
to globally regular gravitating monopoles ||f|,|5],|5| which exist up to a maximal value of the 
gravitational coupling. For higher values of that coupling, the Schwarzschild radius of the 
solution becomes larger than the radius of the monopole core. 

Considering only the theory with a scalar Goldstone field leads to a different type of 
topological defect J7|, the so-called global monopole. Like all global defects this has infinite 
energy resulting from the \ fall-off of the energy density. Coupling to gravity |||9] leads to 
the observation that the effective mass of the system becomes negative. 

Recently, a self-gravitating magnetic monopole in the spacetime of a global monopole has 



been considered HIT||Pl| . In both papers the potential is the sum of the Higgs potential and 
the analog Goldstone field potential. Thus, the interaction between the global Goldstone 
field and the local Higgs field is only indirect, namely through the coupling to gravity. 

Considering the composite topological defect, the effective mass was found to be positive 
or negative, depending on the coupling constants of the model. 

In this paper we continue the investigation of this system, allowing a direct interaction 
between the matter fields. This extra interaction is implemented by adding to the potential 
a gauge invariant term as follows : 

W\x a ) = y(0> a -^)(xV-% 2 ) • (1) 

Here, we are mainly interested in the analysis of the critical behaviour of the composite 
system considering now the most general, gauge invariant potential. Because of the high non- 
linearity of the set of coupled differential equation, an analytical analysis becomes impossible 
and only a numerical analysis can provide the results. 

This paper is organized as follows. In Section II we describe our model and the Ansatz. In 
Section III, we give the equations of motion, the boundary conditions and the analysis of the 
asymptotic behaviour of the Goldstone and Higgs field functions. We present our numerical 
results in Section IV and describe how the extra term (fj) in the potential provides new 
results concerning the behaviour of the fields near the defect's core as well as concerning 
the effective mass of the system. We observe e.g. a strong dependence of the mass on the 
coupling constant A3. We give our conclusions in Section V. 



II. THE EXTENDED MODEL 

This model is described by the following action which is composed of the action for the 
gravitating global monopole and the action of the gravitating local monopole: 
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S = S G + S M = J CcV^gd 4 x + J CmV-9^x , (2) 
with the gravity Lagrangian Co'- 

and G denotes Newton's constant. 

The matter Lagrangian Cm of the extended model with extra direct interaction between 
the global Goldstone field x a an d the Higgs field <ft a reads (a = 1, 2, 3): 

£m = -\F; v F^ a - l -{D,<p a ){D^ a ) - \(d,X a Wx a ) ~ V(<f> a , X a ) , (4) 
with covariant derivative of the Higgs field 

D^ a = d^ a ~ ee abc Al<f) c , (5) 

field strength tensor 

F% = d,A a v - d v A* - ee ahc A\A c v , (6) 
and e being the gauge coupling constant. The potential V(<p a , x a ) is given by: 

v(<P a , x a ) = y - vlf + j (xV - vlf 

where the third term on the rhs couples the two sectors directly to each other with coupling 
constant A3. Ai, A2 denote the self-coupling constants of the Higgs and Goldstone field, 
respectively, while rjx, r/ 2 are the corresponding vacuum expectation values. 

The potential (0) has different properties according to the sign of A = AiA 2 — A3. For 
A > the potential has positive values and its minima are attained for <^ a — 77^, Xa — V21 
for which ([?]) is obviously zero. For A < 0, these configurations become saddle points and 
two minima occur for 

0^ = 0, xI = vI + yv1 (8) 

^3 



and 



Xl = 0, ti = V 2 2 + Y^- (9) 
^3 



The potential's values for these extrema are, respectively, 



4 4 

V ■ = —A V ■ = —A (10) 

' mm a\ ' mm , \ ^— 1 V w / 

4Ai 4A2 



which are negative since A < 0. 
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A. The Ansatz 

The Ansatz for the metric tensor in Schwarzschild-like coordinates reads: 

ds 2 = -A 2 (r)N(r)dt 2 + N-\r)dr 2 + r 2 (d9 2 + sm6 2 dtf) (11) 
where we define for later convenience the mass function m(r) as follows : 

N(r) = 1 - 2a 2 q 2 - . (12) 



The Ansatz for the global Goldstone field x a , the Higgs field <p a and the gauge field in 
Cartesian coordinates reads: 

<p a (x) = V ih(r)x a , (13) 

X a (x) = Vl f(r)x a , (14) 

er 

and 

A a (x) = . (16) 
Substituting the above configurations into the matter Lagrangian density we obtain: 

C M = -Ati dr r 2 A[NlC(f,h,u)+U(f,h,u)} , (17) 
Jo 

where 



and 



+ ^C> 2 " 1)(/ 2 " f) ■ (19) 

The prime denotes the derivative with respect to r. 
The gravity Lagrangian C G is given by 

1 r°° 

Cg= 2GJo drr ( N -^ A '- ( 2 °) 
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III. EQUATIONS OF MOTION 



Varying (0) with respect to the matter fields and gravitational fields and introducing the 
dimensionless variable x and dimensionless mass function fi(x) : 



x = erjir , fi(x) = erjim(r) 
we obtain the following set of differential equations: 



d 

dx 



x z AN 



dx 



A 



2f + x 2 Plf(f-q 2 ) + x 2 p!(h 2 -l)f 



(21) 



(22) 



d_ 

dx 



2 a * T dh 


= A 


x 2 AN— 


dx 





d 

dx 



AN— 

dx 



2u 2 h + x 2 f3l(h 2 - l)h + x 2 plh(f 2 - q 2 ) 
u(u 2 - 1) 



A 



+ uh 2 



x- 



(23) 



(24) 



with 



and 



with 



dx 



[xAN) = [1 - 2a 2 x 2 U]A 



U 



(u 2 -l) 2 u 2 h 2 f 2 f3 2 



2x 4 



■f ^r- + ~, + ^{h 2 - I) 2 + ^{f - q 2 ) 2 



x 2 x 2 ' I 



+ ^{h 2 -l){f-q 2 ) 



dA 9 , - 
— = 2o?AxK 
dx 



2 \ dx I 2 \ dx j x 2 \dx I 



The equations only depend on the dimensionless coupling constants: 

a 2 = AnGrj 2 , (3 2 = A^/e 2 , i = 1, 2, 3 , q = r? 2 A?i • 
With the definition fll2"|), Eq. (^) can be brought to the form 



fi' = a 2 x 2 (K+(U-^-)) . 



x A 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



The finite energy of the solution can then be obtained by taking the value of fi(x) at infinity. 
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A. Boundary conditions 

The boundary conditions at the origin which follow from the requirement of regularity 
read: 

u(i = 0) = l, f{x = 0) = , h(x = 0)=0 , fi(x = 0) = 0. (31) 

In fact the behaviour of the function fi(x) near the origin is p(x) ~ —a 2 q 2 x. The requirement 
of finite energy solutions leads to: 

f(x = oo) = q , h(x = oo) = 1 , u(x = oo) = , N(x = oo) = 1 — 2a 2 q 2 . (32) 



B. Asymptotic behaviour 

The integration of the equations for generic values of the parameters needs a better 
understanding of the asymptotic behavior of the functions f(x) and h(x). Two different 
types of behaviour for the functions / and h in flat space seem possible. Either: 



or 



/(») = , + $ + £>(*) , B^ m j^- W) (33) 



h(x)= 1 + Ciexp(pix) + C 2 exp(-pix) + C 3 exp(p 2 x) + C 4 exp(-p 2 ^), 

f{x)= q + Cx exp(pix) + C 2 exp(-p x x) + C 3 exp(p 2 x) + (7 4 exp(-p 2 x) , (34) 

where p 2 , p\ are the eigenvalues of the matrix 

(qk Q'k) ' (35) 
Our numerical analysis strongly suggests the following results : 

• for PfP% — /3| > the solutions obey the asymptotic behaviour (0) for / and h, 

• for — /3f < the functions / and h have the asymptotic behaviour (|34]) . In this 
case, however, one of the eigenvalues of (^) becomes negative, consequently its square 
root is complex and the functions /, h thus oscillate. This behaviour is clearly observed 
when we set fl\ = 0, j3 2 7^ and increase (3$ from 0. For @ 3 = 0, the solutions exist 
and /, h increase monotonically, as soon as (3$ ^ 0, however, oscillations occur. 
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IV. NUMERICAL RESULTS 



Because of the reasons given previously, we restrict our analysis in the following to the 

case PlPl - PI > 0. 



The limit /3 3 = was studied in detail in [ 10 1 , [|TTJ. Here we discuss how the new term 



influences the behaviour of the solutions. One of the main features is that the Higgs function 
h(x) does not reach its asymptotic value h(x = oo) = 1 monotonically. It first reaches a 
maximum h max > 1 for x < oo and then decreases to 1. This phenomenon, which can be 
expected from the inspection of ( |3~3"D since A > for /3 3 ^ 0, is illustrated in Fig. 1. This 
is different from the phenomena observed in [[HJ, [ IT| . There, for all values of the coupling 



constants, the function h(x) was observed to be monotonically increasing from to 1 as 
indicated in Fig. 1 for /3 3 — 0. The Goldstone field functions f(x) reaches its asymptotic 
value q for increasing values of the coordinate x when /3 3 is increasing. This again can be 
explained by (^) since the value B is a decreasing function of /3 3 with a sharp drop at 
j3 3 ~ 1. Thus, the function f(x) decays less strong for higher values of /3 3 . 

At the same time, the minimum of the metric function N decreases with increasing /3 3 , 
while the matter function fi reaches a maximum at roughly the same x at which the function 
N attains its minimum and than drops down to its asymptotic value which determines 
the mass of the solution. Clearly, for (3 3 = 0, this asymptotic value is positive, while for 
increasing /3 3 , it decreases and becomes negative for large enough j3 3 . 

Another feature of the new term is that, for all parameters but /3 3 fixed, the classical mass 
of the solution decreases when /3 3 increases. This is illustrated by means of Fig. 2 where we 
have plotted the evolution of the mass as a function of /? 3 for three different combinations of 
the coupling constants a, q. For q = 1.0, a = 0.4, the mass of the solution is already negative 
for (3% = indicating that the influence of the global monopole is already dominating in the 
limit of vanishing direct interaction of the Higgs field and Goldstone field. For smaller q, the 
mass becomes negative at some finite value /3 3 = /3 3 . For fixed q, this value is increasing for 
decreasing a. Since the local and global monopole are still "indirectly" coupled over gravity, 
a stronger gravitational coupling, of course, couples the two objects in a stronger way. For 
small gravitational coupling, (3 3 thus has to be raised further to make the influence of the 
global monopole dominating. When the combination f3f/3% — /3| becomes negative, the mass 
of the solution reaches — oo, independent of the combination of q and a. This again can be 
related to the fact that we observe the solutions to become oscillating for \3\fi\ — /3| < 0. 

We also studied the way the solution bifurcates into a black hole when the parameter a 
increases while the others are fixed. Fixing q = 0.4 we have analyzed the critical behaviour 
of the solution and checked that, like for the case /5 3 = 0, the solution bifurcates into a black 
hole for a finite value of a, say a = a c . As demonstrated in Fig. 3 for f3 1 — f3 2 — 1, Q — 0.4 
and /3 3 = 0.8, the function N develops a minimum which becomes deeper while a increases 
and becomes zero for a = a c . The limiting solution thus represents an extremal black hole 
solution with horizon Xh- In Fig. 3, we also show the evolution of the mass with a. Finally, 
the critical solution corresponding to (3 3 = 0.8 and a c ~ 1.085 is displayed in Fig. 4. This 
figure clearly suggests that the limiting solution is not an abelian black hole for x > Xh 
like e.g. in the case of the pure local monopole f|, where the solutions bifurcate with the 
branch of Reissner-Nordstrom (RN) solutions and consequently the functions reach their 
RN values for x > x^. Here, the function / is equal to zero for < x < Xh and non-trivial 



7 



for Xh < x < xq. Similarly, the functions A and h are non-trivial for x < xq and are equal to 
their asymptotic values for x > x , while the gauge field function u reaches its asymptotic 
value for x ~ x^. This solution thus represents a "black hole inside a global monopole" as 
was observed previously for the (3 3 = limit |TT| . 



V. CONCLUSIONS 

In this paper we have analyzed the composite system of a global and a local gravitating 
monopole considering the most general gauge-invariant potential. This potential contains a 
direct interaction between the Golstone and the Higgs field. This term leads to important 
consequences concerning the local behaviour of the fields themselves as well as concerning 
the global properties of the system. One of the most relevant consequences is related to 
the effective mass associated with the composite topological defect. The numerical results 
show a strong dependence of this mass on the coupling constant A3. Although the Goldstone 
and Higgs fields are indirectly coupled through gravity, the extra direct interaction is more 
effective. Increasing the parameter A3, the mass becomes negative, indicating the dominance 



of the global sector over the local one. Compared to the results of the A3 = case pO| . pT 
we observe the modulus of the negative mass to become very large in our system. Another 
point which deserves to be mentioned is that the extra direct interaction term is not positive 
definite. Denoting by Xh=i the value of x for which the Higgs field function h is equal to 
one, h becomes bigger than one for x > Xh=i and the new term in the potential becomes 
negative. However, for a particular choice of the self-coupling constants such that they fulfill 
A = A1A2 — A| > 0, the total potential is positive, vanishing only at the minima (j) 2 a = rjf, 

xl = vl 

As possible extensions of the model studied here, let us mention the coupling to a scalar 
dilaton which arises naturally in low energy effective actions of string theory. The gravitating 



local monopole was studied recently coupled to a dilaton and it was found that in the 
limit of critical gravitational coupling, the solutions bifurcate with the branch of extremal 
Einstein-Maxwell-dilaton (EMD) solutions which are associated with naked singularities. It 
would be interesting to see what sort of critical solution the composite system of a global 
and local monopole reaches since our analysis indicates that the behaviour of the functions 
close to the core of the local monopole is strongly influenced by the global monopole. 



S 



REFERENCES 



[1] P. A. M. Dirac, Proc. Roy. Soc. (London) A 133, 60 (1931). 
[2] G. 't Hooft, Nucl. Phys. B79 (1974), 276; 

A. M. Polyakov, JETP Lett. 20 (1974), 194. 
[3] P. van Nieuwenhuizen, D. Wilkinson and M. Perry, Phys. Rev. 13 (1976), 778. 
[4] K. Lee, V.P Nair and E.J. Weinberg, Phys. Rev. D45 (1992) 2751. 
[5] P. Breitenlohner, P. Forgacs and D. Maison, Nucl. Phys. B383 (1992) 357; 

P. Breitenlohner, P. Forgacs and D. Maison, Nucl. Phys. B442 (1995) 126. 
[6] P. C. Aichelburg and P. Bizon, Phys. Rev. D48 (1993), 607. 

[7] A. Vilenkin and E. P. S. Shellard, Cosmic Strings and Other topological defects, Cambridge 
University Press, 1994. 

[8] M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63 (1989), 341. 

[9] D. Harari and C. Lousto Phys. Rev. D 42 (1990), 2626. 
[10] J. Spinelly, U. de Freitas and E. R. Bezerra de Mello, Phys. Rev. D66 (2002), 024018. 
[11] Y. Brihaye and B. Hartmann, Phys. Rev. D66 (2002), 064018. 
[12] Y. Brihaye, B. Hartmann and J. Kunz, Phys. Rev. D65 (2002), 024019. 



9 




FIG. 1. The metric functions N, /i, the Higgs field function h and the Goldstone field function 
/ are shown as functions of the dimensionless variable x for g = 0.4, a = 0.6, j3\ =02 = 1 and 
three different values of (3%. 
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FIG. 2. The mass in units of -^p- is shown as function of (3$ for three different combinations 
of the coupling constants a and q and for (3\ = 02 = 1- 
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FIG. 3. The mass in units of ^ and the minimum N m of the metric function N(x) are shown 
as function of a for (3% = 0, respectively /?3 = 0.8, and q = 0.4, 0i = 02 = 1- 
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FIG. 4. The profiles of the functions N,A,u,f,h are shown for a « a c = 1.08, /?i = /?2 = 1, 
/% = 0.8 and q = 0.4. 
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